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Abstract 

By a reformulation of the loop expansion in the Resonance-Spectrum-Expansion ampli- 
tude for meson-meson scattering, in terms of s-channel exchange of families of propagator 
modes, we obtain a formalism which allows for a wider range of applications. The con- 
nection with the unitarized amplitudes employed in some chiral theories is discussed. We 
also define an alternative for the Regge spectra and indicate how this may be observed in 
experiment. 

1 Introduction 

In Ref. [1] a nonrelativistic Schrodinger model was proposed for the scattering amplitude of 
non-exotic multi-channel meson-meson scattering, which allows an exact solution in the form of 
an analytic expression for the scattering amplitude. Bound states and resonances are obtained 
through the coupling of the two- meson system to a harmonic oscillator (HO), the oscillator 
frequency being independent of flavor. By flne-tuning the intensity of the coupling, one can 
transform the oscillator spectrum into the spectrum of mesons for all possible flavor combinations. 
The model's results, in particular concerning the resonance and bound-state pole structure in the 
scalar-meson sector [2,3], are well known. Here, it is our aim to show that this model corresponds 
to s-channel exchange of families of propagator modes, similarly to the exchange of a family of 
leading and daughter Regge trajectories. 

In the naivest quark-model picture, quarks and antiquarks are assumed to be conflned to 
a small region in space by strong forces. The bulk of the interactions of the quarks with the 
glue is contained in effective, or constituent, quark masses, whereas the remaining dynamics is 
described by a conflning potential of sorts [4,5]. The quantum numbers of the effective qq system 
are given by the total qq spin S, the relative qq orbital angular momentum L, and the qq radial 
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excitation A^. For the lowest radial states of different flavor combinations, i.e., having = 0, 
the confinement-model parameters can be adjusted to experiment so as to obtain reasonable 
results [6]. However, for higher radial excitations, the results are poor [7]. Of course, the reason 
is the absence of meson loops, in combination with a wrong fine-tuning of the model parameters 
of the = states to the lowest states in the experimental spectrum [1]. 

More elaborate quark models do include meson loops [1,8-10], and predict physical resonances 
instead of fictitious real meson masses only. On the other hand, in an effective, non-microscopic 
picture. Oiler and Oset [11] dynamically generated the low-lying scalar mesons /o(600) (alias a 
meson), /o(980), and ao(980) (the ^0(800) or n meson was not studied), by means of coupled- 
channel Lippmann-Schwinger equations with meson-meson potentials resulting from the lowest- 
order chiral Lagrangian. Their procedure boils down to using an effective four-meson vertex, 
summing up the bubbles from the meson loops, and unitarizing the resulting scattering amplitude. 
Hence, they nicely showed that it only takes a four-point interaction to generate scattering 
poles associated with light scalar mesons. Earlier, Tornqvist had made a similar suggestion [9], 
but went one step further by also assuming pre-existing mesons, which correspond to an input 
spectrum for hare mesons. Unfortunately, his proposal ignored the issue of Adler zeros, and so 
did not allow to conclude that in his model actually two scalar-meson nonets can be generated. 
Later, in Ref. [12] he did include Adler zeros, but in the subsequent analysis, together with 
Roos [13], the Kq(800) (alias k meson) was still not found, probably due to the use of an 
unphysical, negative Adler zero in the 7 = 1/2 case [14]. The latter pole of the isodoublet 5- wave 
scattering amplitude was then indeed generated by Oiler and Oset, together with Pelaez [15], in 
an alternative coupled- channel unitarization scheme. Furthermore, in Ref. [16] Oiler and Oset 
reported the observation of a second nonet, associated with pre-existing meson states, in an 
N/D approach to unitarization. This issue was thoroughly studied by Boglione & Pennington 
in Ref. [17], who came to the conclusion that it is indeed possible to generate two states starting 
from one bare seed (or pre-existing state) only, and that it even may be plausible to dynamically 
generate many states with the same quantum numbers but different masses. That is exactly 
what had been proposed almost two decades earlier in Ref. [2]. 

In Ref. [1] an infinity of seeds was introduced, all related through one parameter, viz. the 
oscillator frequency cj, being the same for all flavors. When, at low energies, the sum is reduced 
to an effective constant, one obtains a four- meson vertex which dynamically generates exactly 
one pole in the case of the scalar mesons. Also, at slightly higher energies, the sum can be 
approximated by its leading term and an effective constant for the remaining sum [18]. Then 
one obtains, apart from the dynamically generated resonance, a second pole associated with the 
one leading seed. However, in general, there is no need to approximate the sum, thus allowing to 
generate several dynamical resonance poles, [2,3,19] and moreover an infinite number of normal 
resonances associated with the infinity of seeds. Furthermore, the model of Ref. [1] can be applied 
to different fiavos, including charm and bottom, with just one set of parameters. The low-lying 
scalar- meson nonet pops up without even being anticipated [2], using the parameters that were 
fine-tuned for the vector and pseudoscalar spectra in Ref. [20]. However, there is an important 
difference with the technique of introducing by hand one or more seeds, namely that the coupling 
constants can be controlled. 

The problem of couplings and propagator modes is not new. In Ref. [21], Delbourgo, Rashid, 
Salam, and Strathdee remarked that it is well-known Hhat Regge trajectories arise from sums of 
infinite sequences of Feynman diagrams in conventional field theory", and, moreover, that this 
"poses the problem of suitable coupling constants". The latter problem was solved in Ref. [21], 
where a strategy was developed for the decomposition of a Regge trajectory in its propagator 
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modes. Here, we shall follow a different but comparable approach. 

In Refs. [1,20], the couplings of the propagator modes were controlled by Clebsch-Gordans 
and 9-j symbols, assuming the ^Pq mechanism [22-25] for quark-pair creation. But this procedure 
led to inconsistencies in the light-meson sector. However, not only the couplings, but also the 
choice of the spectrum of propagator modes is of importance. In Refs. [1,20], such a choice 
was made involving one free parameter, viz. the universal level splitting, which appeared to be 
largely flavor independent, as we shall further elaborate in Sect. [3] below. Based on the choice 
of the spectrum, the couplings issue was solved by the ^Pq recoupling strategy of Refs. [26-28] , 
leaving only one parameter free, which represents the probability of quark-pair creation. As a 
consequence, there are only two free parameters, no matter if the whole sum of propagator modes 
is approximated by just one point vertex, or by the sum of a few propagator modes and one point 
vertex. 

The Schrodinger equation for the meson-meson scattering model was solved with the wave- 
function approach (WFA) in Refs. [1,20]. Later, in Ref. [29], the same dynamical equations were 
solved by an iterative method, which we refer to by Resonance-Spectrum Expansion (RSE). The 
latter method allows more easily a comparison with current theories for meson-meson interactions. 
However, it should be kept in mind that both approaches, WFA and RSE, lead from the same 
dynamical equations to the same expression for the meson-meson scattering amplitude. 

In Refs. [30,31], besides the usual four-meson interaction, a fictitious bare scalar k meson 
was introduced, in a similar fashion as suggested by Tornqvist [9] and Oiler & Oset [16], leading 
to a two-pole description of the Kn 5- wave interaction for energies up to about 1.6 GeV. The 
behavior of these poles, as a function of one overal coupling parameter, is comparable to the pole 
movements of the two lowest lying k (or Kq) poles described in Refs. [2,18]. Namely, the lower 
pole, which is dynamically generated by the Kti-Ktt vertex and associated with the i^'Q(800), 
moves towards larger negative imaginary energies, away from the real axis, according as the 
coupling is decreased, whereas the higher pole, stemming from the seed, approaches the real axis 
for decreasing coupling, ending up at the mass of the fictitious meson for vanishing coupling. This 
similarity in pole behavior inspires us to revisit the meson- meson scattering model of Ref. [1], 
and relate it to the exchange of Regge propagators [21]. 

2 The scattering amplitude 

We define the amputated amplitude for non-exotic two-meson scattering by 

t = v + vQy + vQyvLV + ... = v [i-VLVY^ , (i) 

where V stands for the RSE propagator and VL for the two-meson loop function. In Ref. [29], it 
was shown that the one-channel RSE expression for the £-th partial-wave two-meson scattering 
amplitude follows from 

Vtip) = — Yl T?i jp — ^^^P") = -2i/iproj£ (pro) /if ^ (pro) , (2) 

where p is the center-of-mass (CM) linear momentum, E{j)) the total invariant two-meson mass, 
ji and h^P the spherical Bessel function and Hankel function of the first kind, respectively, the 
reduced two-meson mass, and ro a parameter with dimension mass~^, which can be interpreted as 
the average string-breaking distance. The coupling constants A and qnl are discussed in Sect. HJ 
where in Eq. ( fT6l) the relation between £ and L = L{i) is expressed. 
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The kernel V, which is graphically represented by Fig. [T], is the result of s-channel exchange 
of a system with internal structure, characterized by resonance modes and their masses E^l. In 
the back of our mind we have, of course, the picture of a resonating quark-antiquark system. But 
this is not of much importance yet. So we consider the exchange of propagators with structure, 
rather than of single and pointlike massive objects, in the very same spirit as a Regge propagator. 

At the vertices, each mode couples independently to the meson pairs. The constants gj\fL 




Figure 1: Born diagram for meson+meson^ meson+meson. The wiggly line represents the 
s-channel exchange of a system with internal structure. 

indicate the intensities of the couplings of the propagator modes to the meson pairs. Now, since 
we consider an infinty of propagator modes, we would end up with an infinite number of coupling 
parameters, with which one could describe any two-meson system. However, this freedom is 
strongly restricted by the ^Pq recoupling scheme [27,28], described in Sect. HI In fact, within the 
^Pq scheme, all but one of the parameters qnl are fully determined. The only remaining free 
parameter is A, i.e., the overall three-meson- vertex coupling constant. 

We thus assume that the RSE model can be formulated in terms of s-channel exchanges of an 
infinite spectrum of propagator modes. The expression for the partial- wave loop function Q^lp) in 
Eq. ([2]) stems from the nonrelativistic two-meson loops of expansion ([1]). We have depicted ^le{p) 
in Fig. [21 In the present formulation we are free to adopt different schemes for the determination 
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Figure 2: Loop diagram for meson+meson^ meson+meson. The wiggly line represents the s- 
channel exchange of a system with internal structure. The mesons within the loop are stripped 
off their internal dynamics and just carry the physical masses and quantum numbers. 

of loop functions. This way, with not too much effort, the formalism can also be extended to 
non-exotic baryon-meson and baryon-antibaryon scattering and production processes. 

Besides the coupling constants, we also assume vertex functions which regularize the loop 
integrations. In the RSE, one employs a local spherical delta-shell with radius tq. In the CM of 
the two-meson system, its Fourier transform is given by a spherical Bessel function, i.e., 

je {kro) , (3) 
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where k stands for the CM loop momentum. In the CM system, we furthermore adopt for 
the two-meson loops the techniques developed by Logunov & Tavkhelidze and Blankenbecler & 
Sugar [32,33], and also follow the covariant prescription formulated by Cooper and Jennings [34]. 
This reduces the four- dimensional integration to a three-dimensional one in the CM frame, which, 
assuming spherical symmetry, takes the form 
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2,,rll k^dk f^^ll'l. , (4) 
Jq p'^ — k^ + te 

where k and p represent relativistic three-momentum moduli, and is given by the on-shell 
relative two-meson momentum in the CM frame. This expression straightforwardly leads to the 
loop function fle{p) in Eq. ([2]). 

One obtains a two-meson scattering amplitude T satisfying unitarity, from the amputated 
amplitude t defined in Eq. ([1]), by 

T = XhX = XV [1 - nV]~'^ X , where X^X = XX^ = ^m{n) . (5) 

Consequently, in the loop-expansion approximation of the scattering amplitude, it is sufficient to 
determine the amputated propagator V and the amputated loop function Q. The latter, which 
is depicted in Fig. [21 consists here of two vertices and two free meson propagators. 

The partial- wave amplitude Te{p) is related to the amputated partial- wave amplitude ti{p) of 
Eq. m by 

Te = XjuXe = X\Vi [1 - VLfViY^ . where X|X^ = X^X] = , (6) 

ensuring unitarity for each partial wave. 

One ends up with an expression for the amplitude of non-exotic two-meson scattering, even in 
the case of many coupled channels, which contains only three free parameters, viz. A, tq, and cj, 
besides the constituent quark masses. The resulting amplitudes, for different flavors and a variety 
of orbital quantum numbers, have been confronted with experiment in numerous publications. 



3 The propagator modes 

Soon after T. Regge noticed, in his famous work [35], that bound states and resonances of the 
scattering amplitude are related to their poles in the complex orbital-angular-momentum plane. 
Chew and Frautschi made the observation that the squares of the masses of baryonic and mesonic 
resonances come out on almost linear Regge trajectories [36]. For mesons this has been explored 
in many models, obtaining relations of the form [37] 

El^ = C^ + 2Cuj(2N + L+^^ , (7) 

for meson masses Ej^^. Here, N and L represent the radial and orbital- angular- momentum 
quantum numbers, respectively, of the qq systems. Furthermore, C is a constant which depends 
on the qq constituent flavor masses, and u; is a universal frequency. However, in Ref. [1] it was 
assumed that, at least for low energies, the trajectories of propagator modes for cc and hh systems 
are linear, not quadratic, in their masses. In Ref. [20], the linear mass dependence was extended 
to the light quarks m, c?, and s, too, while in Ref. [2] the linear propagator modes were shown 
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to explain the scalar resonances in non-exotic S'-wave meson-meson scattering as well as the 
scattering data for energies up to about 2 GeV in the Kn + Kr] + Kr]' complex [19,38]. 
For low excitations, the linear mass relation is obvious from Eq. ([7]) if 

C:^u , (8) 

in which case one finds 

Eml-C + u(^2N + L + ^^ , (9) 

but less obvious when C and u are of the same order of magnitude. Now, for u one finds 
in Ref. [20] the value 0.19 GeV, whereas for bb and cc the values 9.45 GeV and 3.12 GeV, 
respectively, are quoted for C, which satisfies well condition ([8]). For the light quarks, C ~ 1 
GeV, which still is not in conflict with condition ([8]) for low energies. Consequently, from the 
good results of the model in Refs. [1,2,19,20,38], one cannot exclude that, at high energies, 
the propagator trajectories will be quadratic in mass. However, since also the results for the 
Kn + Kr] + Krj' complex support linear relations, we will stick here to relation (Q for the 
model's Regge trajectories. 

In previous work on the RSE, we referred to the modes of the exchange propagator ([2]) 
as the confinement bound states. This picture has not been completely abandoned here. On 
the contrary, it will play an important role in determining the parameters that describe the 
propagator modes. Mode masses depend in the first place on the quark flavors flowing in the 
propagator. That information is contained in the constant C of Eq. (Q. However, we will argue 
in the following that the mode level splittings are largely flavor independent. 

The mesonic resonances extracted from experiment are organized by flavor content, J^'-'I'^ 
quantum numbers, mass and width. Based on the bb and cc spectra, it was concluded in Ref. [1] 
that, in principle, there must exist an infinity of such states, though most of the excited states are 
difficult to observe because of the many open two-mesons channels to which they couple, albeit at 
higher masses obscured from observation because of the many two-meson systems which couple 
to qq. Accordingly, we expect an infinite number of scattering poles in meson-meson scattering, 
here represented by 

E = Po, Pi, P2, ... . (10) 

Unitarity then requires that in the one-channel restriction, assuming the poles fllOl) to be simple, 
the elastic scattering matrix S be given bj|3 [9] 

(E - P*) (E - P*) {E - p;) . . . 

SiE) = . (11) 

(E-Po) iE-Pi){E-P2)... 

If we suppose that the resonances ( fTO!) stem from the spectrum of modes of propagator ([2]), 
given by the real quantities 

E = Eq , El , E2, ... , (12) 

then we may represent the differences (P„ — E^), for n = 0, 1, 2, . . ., by the complex mass shifts 
AEn- Thus, we obtain for the unitary S-matrix the expression 

{E-Eo- AEo*) {E-Ei- AE^*) {E-E2- AE2*) . . . 

S{E) = . (13) 

{E-Eo- AEo) {E-Ei- AE^) (E - E^ - AE^) . . . 



Note that we do not consider here a possible overall phase factor representing a background. 
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So we assume here that resonances occur in scattering because the two-meson system couples 
to certain modes of the propagator ([2]), usually of the qq type, viz. in non-exotic meson- meson 
scattering. Let the strength of the coupling be given by A. For vanishing A, we presume that the 
widths and real shifts of the resonances also vanish. Consequently, the scattering poles end up 
at the positions of the mode spectrum (IT^ . and so 

AE„ — ^ for n = 0, 1, 2, ... . (14) 
AiO 

As a result, the scattering matrix tends to unity, as expected in case there is no interaction. 
The scattering amplitude ([1]) exactly satisfies these requirements for the propagator and loop 
function ([2]). As a consequence, it seems that one may only deduce an approximate mode 
spectrum from experiment. Its precise masses Enl can then be found by comparison to scattering 
and production data, once the full scattering amplitude has been composed. However, in Sect. [7] 
we shall see that in production processes the Regge spectrum may become visible. 

In order to set out with the task to find a reasonable ansatz for the mode spectrum of 
our propagator, let us assume that the spectrum of mesonic quark- ant iquark systems can be 
described by flavor-independent HO conflnement. Then, for each pair of flavors, an inflnite set 
of mesons exists with all possible spin, angular, and radial excitations. But unfortunately, for 
most flavor pairs only a few angular and even fewer radial recurrencies are known [39]. If we 
do not distinguish up and down, but just refer to non-strange (n) quarks, then we have at our 
disposal four different flavors: n, s, c, and h. These can be combined into ten different flavor 
pairs, each of which may come in two different spin states: or 1. This gives rise to, in principle, 
twenty different meson spectra. With some 150 known mesons, this means 7.5 angular plus radial 
excitations on average, per flavor pair. This is much less than e.g. the known excitations of the 
positronium spectrum. No wonder that it requires some imagination to guess economic strategies 
for the description of mesons. 
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Figure 3: Nonstrange, charmonium, and bottomonium J^*-^ = 1 states compared to the corre- 
sponding states from a harmonic-oscillator spectrum. The level spacing for the oscillator equals 
0.38 GeV. 

As one may verify from the latest Review of Particle Physics [39], known vector states are 
more numerous than any other type of mesonic resonances, since they are easier to produce. 
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Consequently, in order to structure the mode masses, we begin with the vector mesons, carrying 
quantum numbers J^*" = 1 . In Fig. ([3]) we compare the observed nn, cc, and bb vector 
states with the HO possible states. Most of the data are taken from Ref. [39]. The p(1250- 
1290) signal was originally reported in Refs. [40-44], and has very recently been confirmed in a 
coupled-channel data analysis [45]. The T{1D) has been observed in Ref. [46]. 

The charmonium vector states, shown in Fig. ([3]), bear many similarities with the two-particle 
HO: a ground state in a cc S'-wave, and higher cc radial excitations that are almost degenerate 
with the cc D-wave states. Also, except for the ground state, the level spacings are roughly 
equal. In Refs. [1,20,47], the mechanism was discussed which turns the HO spectrum into the 
charmonium spectrum, also including the ground-state levels. 

For the p and T vector states, also shown in Fig. ([3]), we see a very similar pattern: the 
qq S-D splittings are slightly larger, while the p(770) ground state of the p spectrum and the 
T(15') ground state of the T spectrum also come out far below the corresponding oscillator 
ground states. From Fig. ([3]) one may moreover conclude that, as far as the level splittings are 
concerned, there is not much reason to separate the light-quark sector from the heavy quarks. 
The mechanism which turns the oscillator states into the p and T resonances is discussed in 
Ref. [20]. 

What we learn from the above comparison is that the level splittings are largely flavor in- 
dependent and more or less constant, whereby the ground-state level of the mode spectrum is 
determined by the effective flavor masses. 
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Table 1: The experimentally observed light positive-parity mesons. 

In Table [H we show the experimental spectrum of light positive-parity mesons. Only for 
the /2 states, almost enough resonance data are available, to allow for comparison with HO 
confinement. We therefore assume that the states in the first two /2 columns contain mostly 
non-strange qq pairs, and in the next two predominantly ss pairs. Furthermore, the quark pair 
may come in a relative P-wave (1*^* and 3^*^ column) or F-wave (2'"'^ and 4*'' column). From the 
values for the central mass positions as given in Ref. [39], we collect in Table [2] mass differences 



8 



for a selected set [48] of /2 states. For the spectra of Fig. Owe deduced a level spacing of 0.38 
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— m 
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0.34 ± 0.02 GeV 



Table 2: The experimentally [39] observed mass differences for isoscalar light positive-parity mesons 
with J = 2. 

GeV, which agrees well with the splittings in Table [2J As a first approximation, it thus seems 
reasonable to adopt for the masses of the propagator modes the expression 

M[fJ-N,L)=mf + mj + u(2N + L + ^ . (15) 

Here, / and / represent the flavors of the quark and the antiquark, mj and mf their respective 
masses, and ui the oscillator frequency. 

In the following, let us study some details of formula (fT5l) The vector-meson states have unit 
total angular momentum, J = 1, and unit qq total spin, S = 1. Hence, since the parity of 
vector- meson states equals P = — 1, their orbital angular momentum can be L = (S-wave) or 
L = 2 (D-wave). From formula f|T5l) we then understand that, for HO confinement, the vector- 
meson states with {N,L = 2) are degenerate with the vector-meson states with (A^ + 1,L = 0), 
as shown in Fig. ([3]). For other flavor and spin excitations similar results emerge. One obtains 
a very regular, equally spaced spectrum of propagator modes, with an oscillator frequency a;, 
which comes out at about 0.19 GeV for the data. 

Non-strange (nn) and strange (ss) configurations double the number of isoscalar states into 
SU (3)-flavor singlets and octets. But one should be be aware that all states are mixed through 
the meson loops. Hence, like in Nature we will not find pure angular, radial, or flavor excitations 
for resonances of the scattering amplitude ([1]). Lattice calculations reveal that it may even be very 
hard to disentangle the various configurations showing up in /q systems [49]. Moreover, meson 
loops influence the precise resonance shapes. Some come out broad, others narrower. Also, the 
central resonance positions may shift substantially (100-300 MeV [1,2,20]) with respect to the 
propagator mode spectrum. 

Now, as discussed in the beginning of this section, we may not exclude the possibility that 
the trajectories for the propagator modes are quadratic in mass. Hence, we may thus very well 
assume that relation f[T51) is only an approximation, valid for low energies, and refer to propagators 
([1]) as Regge Green's functions [21]. We leave the study of a more precise relation with string 
theory [50] for future research. 
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4 Vertices 



In this section we study how the propagator modes couple to the meson pair. Thereto, we 
characterize each meson by its quantum numbers. This is in part guess work, since we only have 
at our disposal the total spin J, the parity P, and, for flavorless mesons, the C-parity. But let 
us suppose here that we also have knowledge of the orbital angular momentum L, the internal 
spin S, and the radial excitation n of the quark pair that constitutes the meson. Hence, a two- 
meson system then consists of the sets of quantum numbers {Ji, Li, Si,ni) and {J2, L2, S2,n2), 
characterizing each meson, and also the quantum numbers describing the relative motion of the 
two mesons, viz. J, i, s, and n. The propagator modes are similarly characterized by a set 
of meson quantum numbers {J, L, S, N). Hence, the complete coupling is given by the matrix 
element of the transition operator O 

(( Ji, Li, ^1, m) ; ( J2, L2, ^2, n2) ; J, s,n\0\ J, L, S, N) . (16) 

Given this form of the coupling constants, it is advantageous to determine the scattering 
amplitude from the partial-wave expansion. In the two-meson CM system, assuming spherical 
symmetry, we define 

T(p) = Y.i2e+i)Peip-p') Up) . (17) 

An important ingredient for the vertices employed in the WFA and the RSE is the ^Pq 
mechanism for quark-pair creation of Micu [22] and Carlitz & Kislinger [23], which has been 
worked out in more detail by Le Yaouanc, Oliver, Pene & Raynal [24], Chaichian & Kogerler 
[25], and later by Ribeiro [26]. A complete expression of the latter matrix elements for all 
possible quantum numbers and different effective quark masses is given in Refs. [27,28]. The 
resulting couplings have been employed for perturbative calculus of particle widths and mass 
shifts in Refs. [51-54]. In the latter works, fiavor mass and the universal oscillator frequency 
were combined to one parameter, different for each flavor. 

In Ref. [29] it was shown that the spectral representation of the Green's function, viz. 



E 



|J-^i(ro)|' 2/i Fi(^,ro)Gi(P,ro) 



E-E^L r2 W{FL{E,ro),GL{E,ro)) ' ^ ' 

connects the WFA expression for the full propagator (righthand side [1]) with the RSE iterative 
result (lefthand side [29]). The full propagator is obtained by constructing an exact solution of 
a nonrelativistic stationary equation for meson-meson scattering containing a confining part He- 
The expansion in eigensolutions of He is obtained by an iterative method. The set of functions 
{^NL ] iV = 0, 1, 2, . . .} represents, for orbital angular momentum L, a full set of radial eigenso- 
lutions, with eigenvalue E^^^, of the confining part He of the full Hamiltonian. Furthermore, Fl 
and Gl represent two linearly independent solutions of H^. for any value of the energy E. 

The mode distribution in the CM frame of the two-meson system is contained in (^o)I^ 
of Eq. (fT8|) . In the RSE, the mode distribution is contained in the matrix elements of Eq. (fT6l) . 
Henceforth, we put aside the WFA description of the propagator modes and its coupling to meson 
pairs, and concentrate on the iterative RSE description. Hence, when we abbreviate the result 
of matrix element (fT6l) for two-meson channel i by gi/{N, L), then we propose here to write for 
the partial-wave Regge propagator connecting the two-meson channels i and j the expression 

4 £ f ^-ANML) ,,,, 



10 



where L = L{i) and in which form the RSE formahsm takes a shape similar to the result of 
Ref. [21]. The sum in L is usually very much restricted, since most couplings gi^e{N,L) vanish, 
as follows from the details [27,28] of Eq. 

The vertices connecting the multi-mode propagator to the two-meson systems, are calculated 
in the CM system. On a basis of HOs, for each mode one assumes a spatial distribution of 
two quarks and two antiquarks, one colorless pair for the propagator mode, and an equally 
colorless ^Pq pair. These spatial distribution functions are decomposed on the basis of two-meson 
distributions [55]. The coefficients of such a decomposition are called recoupling constants [27,28]. 
In Table [3] we have collected some of the ^Pq recoupling constants, merely as a demonstration. 



vertex 

PC ^^jPC^^jPC 



recoupling coefficients 
{g{n)y X 4" 



nO- 



or 
or 

or 



oO- 
or 

ol" 

or 
or 



0,0 
0,0 
0,0 

1,1 

ri 



114(2^^-3)^ 

^(2n-3)2 
^(2n + 3) 
^(2n + 3) 



nO- 

nO 



++ 



oO- 
oO- 
lO' 
oO- 

or 
or 
or 
or 
or 

oO^ 

or 

ol^ 



+ 00-+ 

+ 10-+ 

+ 10-+ 

+ 01++ 
+ 01'" 
+ 01"" 

+ ol^=2 
+ 0I+" 
+ 0O++ 
+ 0I++ 
+ 0I+- 



0,0 
0,0 
0,0 

1,1 

0,0 
2,2 
0,0 
0,0 

1,1 

0,0 
0,0 
0,0 



6x242 



2i(^+l) 

_i_(2n + 3)(n-l)2 

n(2n+ l)(2n + 3)(n-3)^ 

j_ 

12 

^(2n + 5) 
■(2n + 3)(n- 1)2 



1 



4x63 ' 

1 r(2n + 3)(2n-5)2 

j_ 

12 



5x63 ' 



r(2n + 3)(n-3)^ 



432 



_i_(2n + 3)(n-2) 



Table 3: Vertex recoupling constants g{n) for the radial excitations n (n = 0, 1, 2, . . .) of pseudoscalar 
and scalar modes of the propagator ([ID , for the case of equal effective quark masses [27, 28] . We have 
characterized the two-meson systems by the quantum numbers r-J^'^' {r for radial excitation), the relative 
two- meson angular momentum i, and the total two- meson spin s. In one case, we have also indicated 
the q-q internal angular momentum of the meson. For the other cases, the lowest possible quantum 
numbers are assumed. 

For n = 0, the recoupling constants squared in one column add up to 1, with the proviso 
that meson pairs with two different mesons count twice, as actually we should have repeated the 
corresponding lines in the table for the interchanged pair. This result reflects the fact that we 
consider a properly normalized distribution for each mode, and an orthonormal set of two-meson 
distributions. 
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For n > 0, the recoupling constants squared in one column do not add up to 1. The reason 
is that for n > more possible two-meson systems couple to the propagator modes. If those 
were included in the tables, we would obtain unity for all n. A full table would have infinite 
length. Moreover, even when limited to J < 10 and only for a few radial excitations, one would 
easily end up with a table of hundreds of pages. However, with a fast computer code in Fortran, 
based on the expressions given in Refs. [27,28], and which moreover takes care of the various 
possible isospin combinations and also allows for unequal effective flavor masses, the absence of 
such tables is no limitation. 

5 Comparison to other models 

In their comment [56] on the work of Tornqvist & Roos [13], Isgur & Speth pointed out that, 
since exotic channels do not couple to the propagator of Eq. ([2]), the corresponding scattering 
amplitude vanishes, which is not in agreement with experiment. Furthermore, they argued that 
in the Jiilich model [57] t-channel processes lie at the origin of a hroad dynamically generated 
pole in the / = S'-wave pion-pion scattering amplitude, and not s-channel propagator modes. 
Finally, they remarked that also the ao(980) and /o(980) are largely due to t-channel forces. 

Now, it seems indisputable that with zero couplings their is no scattering. Nevertheless, the 
first observation of Isgur & Speth is incomplete. 

Namely, when just arbitrary couplings and seeds are involved, then there is no clear pre- 
scription how to handle exotic channels. However, when propagator modes and couplings are 
related, like in the choice of Ref . [1] and the corresponding recoupling constants of Ref . [28] , then 
quark-interchange processes, and so scattering in exotic channels, make part of the interactions 
that can be handled in the model. 

In Ref. [28], the quark-interchange diagrams were suppressed in order to single out the effect 
of the propagator modes in s-channel exchange via ^Pq pair creation. Consequently, although in 
principle there is no limitation in the model inhibiting the study of exotic channels, it is just 
neglected in the RSE approach to non-exotic meson-meson scattering. A study also including 
the quark-interchange diagrams has been carried out by Bicudo and Ribeiro in Ref. [58]. 

Moreover, in Ref. [59] Fariborz, Jora & Schechter observed that neglecting the t-channel 
exchange of a p meson in J = S'-wave tttt scattering does not remove the a pole from the 
scattering amplitude of their model, and only shifts it in a modest way. Also note that the 
dynamical "cr(400)" pole generated by the Jiilich group [57] is considerably lighter and broader 
than generally accepted. 

In Ref. [17], also Boglione & Pennington remarked: "s-channel dynamics is not all that 
controls the scattering", when refering to Ref. [13]. However, one may not conclude from spicing 
up with t-channel exchange, a model which only accounts for one (or a few) propagator modes, 
moreover in the ladder approximation, that meson-exchange contributions are really needed. In 
their replies to the comments on their work, Tornqvist and Roos stressed the concept of duality. 
Hence, it is not even clear whether t-channel exchange should be considered at all, when all 
propagator modes in s-channel exchange are accounted for. Duality [60, 61] probably just works 
the other way around as well. So either take all t-exchange contributions into account, ideally in 
an untruncated fashion and not just in the ladder approximation, or all s-exchange contributions, 
but not both. 

The behavior of the /o(980) and ao(980) poles, for variations in the RSE couplings, was 
studied in Refs. [62,63]. No further t-channel exchange was needed in order to describe these 
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scalar resonances by the lowest dynamically generated poles of the model in KK, for respectively 
/ = and 1 = 1. The modelling is just a bit more complicated than for the a and the k, because 
lower lying, but not very strongly coupled, channels are involved as well, namely vrvr and ?77r, 
respectively. In fact, the correct way to describe the /o resonances is, of course, by employing 
a coupled- channel approach to the / = S-wave [2,64]. Then one obtains in the RSE both the 
/o(600) (a) and /o(980) resonances, in agreement with data, besides other scattering observables 
like phase shifts, line-shapes, and inelasticities [64]. 

Recently, the true pole positions of broad states were highlighted [17,65,66]. However, pole 
positions will always depend on the properties of a specific model. For narrow resonances there 
will be agreement, to some extent, among different models. But for broad structures it is rather 
unlikely that different appoaches will lead to exactly the same pole positions. Boglione & Pen- 
nington wrote: "fitting data along the real axis cannot accurately determine the true pole position 
of a broad state without an analytic continuation, or a very specific model" [17]. Hence, unless 
we fully agree on the perfect model to describe the available data, true pole positions for broad 
structures do not exist. At best, we may agree on whether a specific pole exists. 

In Ref. [18], for the description of elastic P-wave nK scattering and the K*{892) resonance, 
the propagator of Eq. fllSp was approximated by the first term plus a constant representing the 
remainder of the sum: 

Such a procedure is equivalent to the approaches of Tornqvist [12] and Oset & Oiler [16], as 
discussed in the introduction. In the latter approaches, however, one looses track of the relations 
among all coupling constants, thus needing to introduce arbitrary parameters, like a and (3 in 
Eq. (120!) . while also contact with different flavors and other angular momenta is lost. Nevertheless, 
in the approximation of Eq. (!20|) . one can study sufficiently well the properties of the dynamically 
generated and the lowest qq resonances. 

In the case of elastic P-wave irK scattering, no dynamically generated resonance is found. So, 
when the overall coupling is decreased, the K*{892) pole returns to Eq [18]. But for non-exotic 
elastic S'-wave ttK scattering, the dynamically generated -K'g(800) resonance appears, besides 
the Kq{1430). We may thus conclude that the contact term indeed absorps those terms of the 
Regge propagator which are not accounted for. In the loop sum ([1]) its contribution is negligible, 
provided a complete Regge propagator is exchanged. Hence, this result seems to suggest that 
the contact term is not necessary at all in a microscopic formulation. 

The effect of hadron loops on the spectra of mesons and baryons has been studied by var- 
ious groups, and for a variety of different confinement mechanisms [8,67-69]. For mesons, the 
procedure usually amounts to the inclusion of meson loops in a qq description, or, equivalently, 
the inclusion of quark loops in a model for meson-meson scattering. This results in resonance 
widths, central masses that do not coincide with the pure confinement spectrum, mass shifts of 
bound states, resonance line-shapes that are very different from the usual Breit-Wigner ones, 
threshold effects and cusps. In particular, it should be mentioned that mass shifts are large and 
negative for the ground states of the various flavor configurations [20]. Unquenching the lattice 
is still in its infancy, at least for the light scalars, as we conclude from Ref. [70]. However, its 
effects should not be underestimated. Hence, ground-state levels of quenched approximations for 
qq configurations in relative S'-waves must be expected to come out far above the experimental 
masses. 
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6 Quark- interchange contributions 



As mentioned in the previous section, the model is not hmited to the study of non-exotic channels. 
However, ^Pq pair creation/annihilation, through which process meson pairs couple to the Regge 
propagator, does not work for exotic channels. Nevertheless, the alternative, which is quark 
interchange, does give contributions to all possible hadronic final-state interactions, hence also 
for exotic two-meson channels. 

Quark-interchange contributions to meson-meson interactions can be determined [58] with 
the very same techniques that were developed in Refs. [26-28]. 

7 Experimental results for the Regge spectrum 

As may be concluded from Eq. ([1]), and from the expressions for Vi and given in Eq. ([2]), 
the dressed partial-wave RSE propagator for strong interactions takes the form (restricted to the 
one-channel case and leaving out some parts not essential for our discussion in this section) 

{oo |5'A''l| 1 
1 - zj, (pro) h'l^ (pro) E ^ ^ • (21) 

n=0 — ^NL I 

This propagator has the very intriguing property that it vanishes for E —>■ E^l. Hence, one 
may wonder what happens in a physical process when the propagator does not allow any signal 
to pass. We shall show in the following that this phenomenon can be, and has indeed been, 
observed in experiment, but not in scattering processes. 

The RSE amplitude for strong scattering is given in Eqs. (fT][5l) . One easily verifies that it 
does not vanish in the limit E — * Ej^^. However, for strong production processes, we deduced 
in Ref. [71], following a similar procedure as Roca, Palomar, Oset, and Chiang in Ref. [72], a 
relation between the production amplitude P and the scattering amplitude T, reading 

Pi = 3iipro)+tTehP (pro) , (22) 

which, using Eqs. f ET]) and ([T]), can also be written as 

Pe = jiipro)U,{E) . (23) 

For the latter expression we find, by the use of Eq. (pTj) . that the production amplitude of Eq. (!23|) 
tends to zero when E -Eatl. This effect must be visible in experimental strong production 
cross sections. 

Actually, the primary question here is not so much if a vanishing qq propagator is observable, 
but rather whether the production amplitude always vanishes when E —>■ E^l- In order to answer 
this, we must return to the results of Ref. [71], where we found, for the complete production 
amplitude in the case of multi-channel processes, that Eq. fl23l) represents the leading term, and 
that the remainder is expressed in terms of the inelastic components of the scattering amplitude. 
The latter terms do not vanish in the limit E — ^ -Eatl? as we have discussed above. Hence, the 
production amplitude only vanishes approximately in this limit, in case inelasticity is suppressed. 

However, there are more questions to be responded with respect to the observability of van- 
ishing propagators. Namely, do processes exist where only one partial wave contributes and 
in which processes other than s-channel exchange do not play an important role? Fortunately, 
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the answer to the latter, pertinent, question can be responded afirmatively, because electron- 
positron annihilation into multi-hadron final states takes basically place via one photon, hence 
with J^^ = 1 quantum numbers. Consequently, when the photon materializes into a pair 
of current quarks, which couple via the qq propagator to the final multi-hadron state, we may 
assume that the intermediate propagator carries the quantum numbers of the photon. Moreover, 
alternative processes are suppressed. 

We may thus conclude that, if we want to discover whether the propagator really vanishes at 
E — Enl, then the ideal touchstone is e"'"e~ annihilation into multi-hadron states. But not only 
do we have at our disposal a wealth of experimental results on such processes, there also exist 
predictions for the values of Enl, with L = or L = 2, given by the parameter set of Ref. [20]. 
As an example, for cc one finds in the latter paper Eq^ = 3.409 GeV and u = 0.19 GeV. Using 
Eq. ( fT5l) . we then get for the higher cc confinement states the spectrum Ei^ = Eq^2 = 3.789 GeV, 
^2,0 = ^1,2 = 4.169 GeV, ^3,0 = ^2,2 = 4.549 GeV, ... . 

The latter two levels of the cc confinement spectrum can indeed be clearly observed in ex- 
periment. For example, the non-resonant signal in e"'"e~ — 7r+7r~'?/;(2S') (see Fig. 5 of Ref. [73]) 
is divided into two substructures [74-76], since the full cc propagator fl?I]) . dressed with meson 
loops, vanishes at E^ = 4.55 GeV [20]. In the same set of data, one may observe a lower- 
lying zero at E2 = 4.17 GeV [20], more distinctly visible in the data on e"'"e~ — *• 7r~^7r~J/ip (see 
Fig. 3 of Ref. [77]). The true cc resonances can be found on the slopes of the above-mentioned 
non-resonant structures [78], unfortunately with little statistical significance, if any [79]. 

In the light-quark sector, where low statistics does not result in sufficient accuracy, one needs 
some imagination for the identification of vanishing amplitudes in e"'"e~ annihilation processes 
into multi-hadron states. Thus, from data given in Ref. [80], one may infer that certain four- 
pion amplitudes vanish at, or near, the predicted [20] 1.097 GeV nn ground state, and near 
the predicted 1.477 GeV nn excited state. Some further evidence, also for the predicted higher 
excited states in nn and ss, is found in data from Refs. [81-85]. 

So we indeed observe minima in production processes, which confirm vanishing qq propagators. 
Moreover, the qq confinement spectrum predicted 25 years ago in Ref. [20] seems to agree well with 
experimental observations for vector mesons. Accordingly, we expect vector-meson qq resonances 
associated with each of the Regge states: one ground state, dominantly in a qq S-wave, and two 
resonances for each of the higher excited Regge states, viz. one dominantly in an S-wave, and 
the other mostly in a D-wave (see Fig. [3]). 

8 Conclusions 

We have found that the meson spectrum can be described by s-channel exchange of Regge 
propagators in non-exotic meson-meson interactions, and furthermore that the bound-state and 
resonance spectrum of mesons is richer than, and different from, the underlying Regge spectrum 
[2,19]. For the latter we generally observe very regular and equidistant level spacings, instead 
of quadratic trajectories. We also have established a link between s-channel exchange of Regge 
propagators and unitarized chiral models constructed for the study of resonances in S'-wave 
meson-meson scattering. Moreover, a method is given to relate the coupling constants and 
seed masses for the latter models. Moreover, we have indicated how the Regge spectrum of qq 
propagator modes may be observed in production processes. 

Finally, the use of Regge propagators for meson physics is seen [1,20] not to be restricted to 
light fiavors, but can, without any further effort, be extended to heavy quarkonia with the same 
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set of parameters, i.e., the effective quark masses, the universal frequency 00, the overall vertex 
intensity A, and the average string-breaking distance tq. 
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